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Relations between the Retention Equations Used in Solid-liquid 
Chromatography Applying Mixed Mobile Phases 
Zależności między równaniami retencji stosowanymi w adsorpcyjnej chromatografii cieczowej 
/ z wieloskladnikowa fazą ruchomą 
Sanucumocra MEN pABNENKNMM YNEPXMBANKA, NPMMENAEMKIAM 
B KMIĄKOCTHO-A4COPÓNMOHHOH xpoMaTorpatpMK C MHOTOKOMIIOHEHTHOŃ NOĄBHKHOŃ cpa3oń 


Among various equations proposed to describe concentration 
equilibria in Solid-Liquid Chromatography systems, there are 
two which have appeared to be pretty useful. One of them is 
the Snyder-Soczewiński equation Ca. 2]: 


: 
= 


in @, - const = ln x, __ (1) 


in which * is the distribution coefficient of a solute k“, 
between the surface and the mobile liquid phase, being a bina- 
ry solution, Further, x is the mole fraction of the solvent 
"2" in this mobile phase, while const. is a quantity independ- 
ent upon the composition of the mobile phase. The theoretical 
considerations which led Snyder and Soczewinski to eqn. 1 were 
different P. J, but both of them concern the case when the no- 
„bile phase is composed of two solvents only. 
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Another useful equation was proposed by Ościk [B]: 


in, = * 4 %, a + 2 as ws + const ye (2) 


where ei o is the appropriate distribution coefficient in the 
pure solvent "i", and ns is the surface excess of this solvent, 
` 4 ne = 9; — * $ | (3) 
Above 0, is the adsorption isotherm of the component "i" on the 
support surface. i ' 

These equations were known for a long time, and numerous 
examples of LC systems were reported, which could be described 
by either eqn, 1 or 2, At the same time it was clear that any 
- transition from one to the other equation is possible, on the 
ground of the theories proposed by Snyder; Soczewiński and Ościk, 
It is the purpose of this paper to investigate this problem in 
more detail. Using methods of statistical thermodynanics and a 
simple lattice model of the surface phase, we shall show that 


both these equations are simplified forms of some more general 
equation. . ga = i 


THEÓRY 


When describing the properties of the surface and the mobile 
liquid phase on a molecular level, we shall use a lattice de- 
scription formalism, which is commonly used in the theories of 
bulk solutions and in the theories of adsorption as well. Thus, 
our present model of LC systems may be described shortly as 
follows: 

1. The mobile liquid phase (multicomponent mixture, in gen= 
era!) may be as well ideal as non-ideal mixture, but still a 
‘regular solution. 

2. The molecules of the solvent mixture (mobile wee) and 
of the solute (analyzed substance ) are adsorbed on a twodimen~ 
sional lattice of adsorption sites, one molecule per one site. 

3. All adsorption sites are energetically equivalent, i. Ges 
the support surface is assumed to be homogeneous, 
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4, It is assumed that small samples of solutes are analy - 
zed, so the concentration of solute in both the mobile and the 
surface phases is small, We shall, therefore, neglect the mu- 
tual interactions between solute molecules in both the mobile 
and the surface phases, 

5, The relatively small solvent molecules adsorbed of dif- 
ferent sites, do not interact with each other. At the śame|time 
the relatively big solute molecules may interact with other 
solvent admolecules. 

Now, consider first the canonical, system partition function 
alu e] T) for the case when the surface phase is ideal, 


$ M! rj Ny (M- ZN) 4 
Ahr) ern e - „An M 6) 


kad. K=q 


Similarly as in our previous paper, M is here the total number 
of adsorption sites, N, is the nymber of admolecules of k-th 
kind, and qq is the molecular partition function załóg: as 
follows: 


qx dyk zk 5 aoe ( zł | A (5) 


where Gy Aye Gz-are the vibrational partition functions of 
the vibrations parallel (q,, q y) and normal (a) to the|surface 
+07 is the partition Fand de for internal vibrations, 
and rotational degrees of freédom of the adsorbed molecules, 
Further, € is the so-called “adsorption energy” which is just 
a difference between the energies of a molecule in the adsor- 
bôd and in the bulk phase, assuming that this molecule is in 
its lowest quantum energy states in respect to all its degrees 
of freedom in the adsorbed and in the bulk phase. 

In case of adsorption from solutions the following nee — 


a 


tions are gia moa ; - 7 


ate | jag: ager 


and, consequently: 
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3 Im Qi 
aj Sy 
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above, pf) is the chemical potential of the molecules of i-th 
kind in the surface phase, Using Stirling s approximations for 
InN, !, and considering that: 


pe ay „ 48) 
- .where pe 45 appropriate chemical pot l in the bulk mobile 
phase, we obtain: ; i í - = 
VVV ' 6) 


-Pk Pa = kT ln = 
3 <A eee ( ) a 


Let us consider now the effects ef the mutual interactions 
between adsorbed molecules. We shall take them into account in 
a simplest possible way using the well-known "Fean-Field Ap- 
proximation”. It is assumed in MFA, that there are no correla- 
tions between adsorbed molecules on the surface. In other words, ` 
the molecules are distributed on the surface completely at ran= -~ 
dom. Tius, in a nearest neighbourhood = a certain admolecuie, 
there should be found on an average (z © N/M) + admolecules 

of i-th kind. Here z © denotes the number of the nearest 
neighbours ~ adsorption sites. Let further we) denote the in- 
teraction energy between two molecules of i-th and j-th kind, 
respectively, which are the nearest neighbours. Then, every Ja, 
adsorbed molecule of k-th kind is subject to the action of a 
mean Toros field Ek 


N. $ 
ED . = w) (- 65 (10) 
: isi i M 


Let the set of subscripts k=i,2,.,.,l. , denote the compo = 
nents pf the mobile phase, whereas the set k= (1+4),....n de- 
note the analyzed substances. According to our assumption (8) 
we have: 


— 
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8, = pts cap ary for k- Q+1),....n (11) 
M WZ . 
and, therefore: 
1 j . 
s 6 ; ; 
AE EW a, Boost T Boeedn (12) 
- | deat 


1 

and EP sz” wę? . X 
t 0 

For the purpose of our present work, we shall neglect the 
interactions between solvent molecules on the surface. at the 
same time, we shall assume that solute molecules may interact 
with solvent molecules on the surface. With this approximation, 
our equations are still valid, agate that molecular dad 
functions (take the form: 


(u) Ek 
qm dw ar fr eP (GE) 2 
7 (i) (zag AA: (13) 
e 


* 1-4 
Consider now the relation between the chemical potential pO. 
activity a? „ and the mole fraction x, of the k-th species in 
the equilibrium bulk phase. We shall accept further a non-symmet= 
rical reference system for standard chemical potentials, in 
which we distinguish “solvents” and the "soluted substances", 
For the “solvents" we have: N 


for kai; isk 


ME ve +-kT ln a,” : Kk51, 2. 1 (14) 


where ph’ is the standard chemical potential of the pure 
solvent . In the same reference system, for "soluted" (analy- 
zed) substances we have: 


p = PA + kT ln Xk ; Erd 45) 


where Pk — is the standard chemical potential of the 


analyzed substance „k“ at. its infinite dilution. Obvious- 
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ly, it is not equal- to m standard chemical. poten = 
tial of the pure solute * From eqns. 9 , 11, 14 and 
15, we obtain: 


be) 


where „je ER ne „O R 2 
< i 3 K 8 i 
175 | x d exp oe 125 (17) 
where : ; pee S 8 * 
a z 4 $ 2 pe) for kul, 2. 71 i 
| b) : EPK (18) 
- Pk.: (b) | 
: 5 YE. - F A += for ke (1+3). sean 
Further, D 23 > ; 
: 8 7 x 7 for k=1,2, os "E 
ab) = ae Se | (49) 
Xk i for k= (11) AT) 
and ? ZŁA: 
"lim ż e WSE ia for 1 2 „ 1 (20) > 
5 a Te Bad “7 | 


In the case when the Hx dae pa is a regular solution, we have 


z 


A, we JĄ pó ZZA ** E 


ist 


where 045 154 is the. excess of the Gibbs free energy of 
mixing oP 17 and "j". In case of a binary mobile phase 


* has the PATRON relation to x : A 
` 5 


kT Ing, = 15 — ker ln J. = 2 5 (22) 


* 


Thus, for the components of the mobile phase, we have 
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(b) 


1 (b) 
— J. d. exp (2) 


121 


(23) 


whereas for the EE substances, we have: 


| 
sę „op Lit ay AA 0.2905 X sd 
ko O) 
£x X. qc EP a 
The „ coefficient > is just the ratio of the con 
centrations (mole fractions) of the analyzed substance in Ans 
surface, and the mobile phase, so we have: 


Pee a (ba) An 


(es) 


Let us consider now the form of the distribution coefficient 
T. when, the mobile phase is the pure solvent, 4“ We shall de- 
note it by® 0) From eqn. 24 , in the limit xmi. we obtain: 


eo [ds (es ES] 


(26) 


So, it can easily be deduced that the distribution coefficient 
Yr (f3) in the mixture {x,} of the solvents, may be expressed 
by means of appropriate distribution coefficients in pure sol- 


vents $ k,O * 


e (x)= ff fea KS a [exp A CZJM 


(e) (27) 
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For the case of binery liquid mobile phases, eqn. 27 gn be 
rewritten to the form: 


cS) (e) 
OC gw tj ao T ES m, ta 2) 
g© 


+ „u SE + m SĘ + Sa x,x, 


X, 


where ne) is the surface excess of the component "i" which, 
for our model of adsorbed phase, tekes the form: 


© . ER E A : ; 
AN, = Tr AERO | e 
and where ` 


z” T ar 0 $ 


RESULTS AND DISCUSSION 


Let us assume that: 

a} the standard chemical potentials SA and Pe are 
equal ; 

b) the interactions in the adsorbed phase bótebóń solute 
and solvent molecules can be neglected, i.e., E> = Ek 5-0 


With these simplifications, eqn. 28 can be SZYN to the 
form: 


ln E. In By - lm (x, X, + nE 94 61) 


+ a X. xa 


Let us assume further that the mobile phase is' ideal, i.e., 
71 = T = 1. Then, from eqn. 31, we obtain: 


10 X = const = ln ( + MOCZ 1)) G2) 
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Let us assume at last, that the second component is 2 more 
strongly adsorbed than the component "i", 1. e., that 10 >i. 
Eqn. 32 reduces then to Soczewiński-Snyder eqn. 1 
(3) 
dk 

ln 1 5 In 17 - in * 2 3) 
In Fig. 1A we have shown the-experimental data by Ościk and 
Chojnacka [4] (benzene~cyclohexane as the mobile phase and Alu- 
mina G as adsorbent) plotted according to eqn. 33. One can 
observe there some serious deviations from a linear behaviour 
vs. Inx,, what means that eqn. 33 is rather a poor approxima= 


~ tion.. 


In Fig. 1 we have shown a) behaviour of the more general 
equation 31. The constant A was assumed to be equal to 
x sx) „and its average . e was used in calculation, found 
Wy averaging the following values in Table 1. 
Table 1 


Solute 2 A; = w 7 a 
Quinoline 15.17 
6-Methylquinoline 14.79 
2,6-Dimethylquinoline 12.22 
Acridine 13.00 


Average 


Looking to the Fig. 18 we can see, that also the more 'gen- 
eral eqn. 31 does not lead to a satisfactory agreement with 
the experiment. Thus, we are going now to investigate the 
applicability of Ościk's eqn, 2. 

This equation is obtained from our general eqn. 28, by as- 
suming that the individual adsorption isotherms of solvents 
8,’ s are linear functions of appropriate bulk activites x, L Se 
We should mention at this moment, that in his original |devel-. 
opment Ościk [3] takes also the differences in bulk and sur- 
face molar volumes of solvents into account. However, it does 
not make any serious difference in our particular case of the 
benzene=cyclohexane mixture. The full Ościk s equation has the 
form: 
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BO Be W Be Jit ag 


Fig. ra A = the experimental data by Ościk and YA Ba LA 
plost d toa rding to e a Aub = the function f(x)= ln XK 
5 + X2 ) ~ ai2x1Xp plotted against x2, . 
provided. thet = 13.8; 21 ~ quinoline, 2 = popula wszy 
3-2 „6-dinothyląvinolino, 4 ~ acridine . 


. 


PA Z ; > AGR (e) © 
In X = x, 1 + x, Xo In 4,2), const nO), 2. 5 64) 
> TM eee 
where ve and v, are kokar volumes of the solute and of the 
0 mixture, whereas Ac 2 is the difference between the 
excesses of the free energies of mixing solvents in bulk and 


surface phase, respectively. Since we assume that the surface 


phase is ideal, this term is very close to our term 22 & 1 * · 
Let Fp (xijdenote. the function, 
> w 29 — | 
FA (x)= In X. * ind: - * In ZĘ - 2 2 * (35) 


` According to Odcik’s eqn. 34, Fk? s should be proportinal to 
nie „In Fig. 2A, we have shown 1 s for the investigated chro- 
matography systems, whereas in B = these functions, normalized 
so that their maxima be still at 0.63. 
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05 
9 3 x — 4 


Fig. 2. A = the functions F,'s from eqn.35 plotted against * 2 
B - normalized functiońs; for explanation see Fig. 1 


“We can see that the normalized excess isotherms overlap each 
other pretty well. This is an evidence that the term Const. nf) 
is very important, i.e., that the interactions between solvent 
and solute molecules on the surface play an important role. 

Of course, one should except that the full equation 28, will 
be the most adequate description of solid~liquid chromatography 
systems. In order to demonstrate it we would have to perform 
calculations of best-fit type, with respect to the two para- 


meters; LG ànd Const = (EP - E,9)/kT. However, we are not 


sure if, having six experimental points, makes such a best- 


fit calculation reliable, - 5 

Moreover, we feel that the present approach should be devel- 
oped. further, to take other important physical factors into 
account, For instance: 7 A 

1} interactions between solvent molecules on the surface; 

2) differences in the surface areas occupied by solute and 
solvent molecules; 

3) surface non-uniformity; 

4) multilayer adsorption, 

We are continuing our investigation along these lines. 


` 
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„STRESZCZENIE 


Wyprowadzono nowe równanie 8 zaléznosé współczynni= - 
ka rózdziału, od składu wieloskładnikowej fazy ruchomej, u 
ględniające nieidealność| fazy ruchomej i oddziaływania czaste- 
czek substancji chromatogtafowanych z cząsteczkami rozpusz= 
czalników w fazie powierzchniowej. Przy pewnych założeniach 
upraszczających z równania tego otrzymuje się dobrze znane 
równania Ościka i Snydera =. Soczewinskiego. : 


x 


Peznuó * 8 


BHBEAEHO - HOBOS ypabreaue -„ONUCHBAKĄCE 38BUCHMOCTK KOSŚÓÓUNUCH- 
Ta pacnhpenenenu OT COCTABA MHOTOKOMNOHGHTEOŃ NOĄBUXHOŃ AZK- 
Sto ypaBHeHWe JUWTHBAET. HEKXGANBHOCTB NOĄBUXHOŃ aan u B38UMO— 
nelicrBue MERAY MOIEKRYJAMM pA3ĄENAGMOTO Beuecrza U pACTBODUTE- 
nel B HernoqBuxHiofi dase. [pu HeKOTOPHX NpeANONOKeHUAX. US BTOTO 
YPABHCHUA MOXHO NOXYUKTB XODOWO e e YyPABHGEKA OcbuuKa 
u Cuafinepa-CO4EBHHCKOTO . 2 


